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We simulated the quantum dynamics for magnetic resonance force microscopy 



(MRFM) in the oscillating cantilever-driven adiabatic reversals (OSCAR) technique. 
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. We estimated the frequency shift of the cantilever vibrations and demonstrated that 

' this shift causes the formation of a Schrodinger cat state which has some similarities 

and differences from the conventional MRFM technique which uses cyclic adiabatic 
. reversals of spins. The interaction of the cantilever with the environment is shown 

to quickly destroy the coherence between the two possible cantilever trajectories. 
We have shown that using partial adiabatic reversals, one can produce a significant 

O" 1 increase in the OSCAR signal. 

> '■ 

• i-H , 

^ ■ PACS numbers: 03.65.Ta,03.67.Lx 

I. INTRODUCTION 

The oscillating cantilever-driven adiabatic reversals (OSCAR) technique is, probably, the 
most promising way to achieve single-spin detection using magnetic resonance force mi- 
croscopy (MRFM) |l|, 0]. In the OSCAR technique, the oscillating cantilever produces an 
oscillating ^-component of the effective magnetic field in the system of coordinates con- 
nected to the rotating rf field, which causes the magnetic moment reversals. In turn, the 
cyclic adiabatic reversals of the magnetic moment cause the frequency shift in the cantilever 
vibrations, which is to be measured. 

While the classical dynamics of the OSCAR has been studied in , the quantum 
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FIG. 1: The geometry of the OSCAR MRFM. rh is the magnetic moment of the ferromagnetic 
particle on the cantilever tip; \x and S are the magnetic moment and the spin of an atom (in the 
ground state); B ex t and B\ are the permanent and rotating magnetic fields. 

dynamics has never been investigated. This paper describes the quantum dynamics of the 
spin-cantilever system in the OSCAR technique. In Section 2, we introduce the quantum 
Hamiltonian and the equations of motion for the system. In Section 3, we present a quali- 
tative analysis and estimates of the OSCAR signal. We show that the OSCAR signal can 
be increased by sacrificing the full reversals of the effective field. Namely, partial adiabatic 
reversals can provide a larger signal than the full reversals. We estimate the maximum 
OSCAR signal generated by a single spin. In Section 4, we present numerical simulations 
of quantum dynamics using the Schrodinger equation. In particular, we demonstrate the 
formation of a Schrodinger cat state which simultaneously includes two possible trajectories 
of the cantilever. In Section 5, we present numerical simulations of the decoherence caused 
by the interaction of the cantilever with the thermal environment. 

II. THE HAMILTONIAN AND EQUATIONS OF MOTION 

We describe the cantilever interacting with a single electron spin 1/2 as a harmonic 
oscillator with the fundamental frequency u c . The dimensionless quantum Hamiltonian of 
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this system can be written (in the rotating coordinate systems) as 

1 



H = ^p 2 + z 2 )+eS x -2rizS z , (1) 



with the dimensionless parameters 



p = P/P , z = Z/Z , e = jB 1 /u c , P = h/Zo, (2) 
dB z 



7] = 



/2{hu c k c ) 1/2 , Z = {%uj c /k c ) 1/2 . 



dz 

Here P and Z are operators of the momentum and the z-coordinate of the cantilever, k c is 
the cantilever spring constant, B>\ is the magnitude of the rotating magnetic field, 7 is the 
magnitude of the electron gyromagnetic ratio, B z is the z-component of the magnetic field 
acting on the spin, including the external permanent magnetic field B ext and the magnetic 
field produced by the ferromagnetic particle. (See Fig. 1.) 

nn 

The quantum dynamics of the MRFM was studied earlier in |3J,|4JJ. In these papers the os- 
cillating z-component of the effective magnetic field acting on the spin was considered to be 
generated by frequency modulation of the external rotating field B\. The z-component of the 
effective field, —2r]z, associated with the cantilever vibrations, was negligible in comparison 
with the z-component of the external effective field. Below we consider the opposite (OS- 
CAR) situation: the oscillating z-component of the effective field is produced by cantilever 
vibrations. 

To simulate the Schrodinger dynamics, we use the dimensionless Schrodinger equation 



z— =H% (3) 



where the wave function ^ is 



^/ = u a (z,r)a + U/3(z,t)P, (4) 

The orbital wave functions u a ^{z,r) have been expanded in the eigenfunctions, u n (z), of 
the oscillator Hamiltonian with time- dependent coefficients. The resulting system of differ- 
ential equations for time-dependent coefficients was solved for a finite set of basis functions 
u n {z). Alternatively, we have found eigenvalues and eigenvectors of the time-independent 
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Hamiltonian Ji for the same finite basis, and constructed from them the time dependent 
solution. These two approaches have been shown to give the same results. 

To simulate the decoherence process, we used the simplest master equation (an ohmic 
model in the high-temperature approximation 



dp s , s 
dr 



1 



2Q 
—2irj(z's 



z — z 



dz 2 
{I 



dz'\ 
' dz' 



-2 U 



(z 2 - z' 1 ) 



(5) 



zs) p S)S i - i^(p s ,-s> - P-s,, i- 

where p SyS i = p s<s i(z, z' , r) is the density matrix, s, s' = ±1/2, D = ksT/huJc is the dimen- 
sionless diffusion coefficient, and Q is the quality factor of the cantilever. 

The function p s , s '(z, z', r) was expanded in the basis of the eigenfunctions u n (z)u* m (z') 
with time-dependent coefficients. These coefficients were calculated numerically for a finite 
set of basis functions. 



III. QUALITATIVE ANALYSIS AND ESTIMATES OF THE OSCAR SIGNAL 

Suppose that initially a single spin is parallel to the effective magnetic field B e ff = 
(e, 0, —2r]z). Under the conditions of adiabatic motion, the spin component along the effec- 
tive magnetic field is an integral of motion. Thus, in the process of motion the spin will have 
the same definite direction relative to the effective magnetic field. Consequently, the orbital 
and spin degrees of freedom are not entangled: the wave function of the system will remain 
a tensor product of the orbital and spin parts. As a result, the spin-cantilever correlators 
are equal to zero, and the Heisenberg equations of motion 

z=p, p= -z + 2r]S z , (6) 

£>x = 2t/zS Z} Sy = 2t]zS x sS z , S z = sSy, 

are equivalent to the classical equations of motion. That is why in this case we can use 
classical estimates for the single-spin OSCAR signal. 

We estimate the OSCAR signal in the following way: under the conditions of adiabatic 
motion (and if the spin is parallel to the effective field) we have 

S z /S = TB z eff /B eff , (7) 



5 



where the upper sign corresponds to the ground state (the spin points in the direction 
opposite to the direction of the effective field). It follows from Eq. (7) that 

S z = ± . VZ (8) 

Now, we substitute this expression into the equation for the cantilever coordinate z 

z + z- 2r)S z = 0. (9) 

Neglecting oscillations with twice the cantilever frequency, we have from Eqs. (8) and (9) 

2n 2 z 

z + zt-i= =0, (10) 

/2r] 2 z 2 m + e 2 



where z m is the amplitude of the cantilever vibrations. Assuming that the cantilever 
frequency shift Au c is small in comparison to the unperturbed cantilever frequency, u c , 
(|Au; c | <C u c ) we have the estimate for the OSCAR signal: 

^= T , i . (id 

Note that the frequency shift (11) decreases with increasing of the cantilever amplitude 
z m . The reason for this dependence is the following. The frequency shift is associated with 
the restoring spin force F s oc z. Roughly speaking, the OSCAR signal is determined by the 
ratio F sm /z m , where F sm is the maximum spin force. When z m increases, the value F sm 
approaches its limit r\ (in dimensionless notation). Thus, the ratio F sm /z m decreases. 

As an example, we estimate the single-spin signal for the experimental parameters in 

k c = 0.014N/m, w c /2tt = 21.4kHz, Z m = 28nm, (12) 
dB z 



1.4 x 10 5 T/m, B x = 0.3mT. 

oz 

Using Eq. (2) we obtain 

7] = 2.9 x 1(T 3 , z m = 8.8 x 10 5 , e = 390. (13) 

It follows from Eq. (11) that the relative cantilever frequency shift is \Au c \/u c = 2.3 x 10~ 9 . 

The condition for the adiabatic motion in terms of the dimensionless parameters can be 
written as 

2 V z m < e 2 . (14) 
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To provide the full reversals of the effective field one also requires 

£ < 2r]z m . (15) 

As an example, in the experiment JlJ we have: 2i]z m = 5.2 x 10 3 , e 2 = 1.5 x 10 5 . Thus, both 
inequalities (14) and (15) are satisfied, i.e. we have full adiabatic reversals. 

Next, we discuss a way to increase the OSCAR signal. It follows from Eq. (11) that one 
should sacrifice the requirement for the full reversals (15) and reduce the amplitude of the 
cantilever vibrations z m and the amplitude of the rf field e, in order to increase the frequency 
shift. The minimum possible dimensional cantilever amplitude can be estimated as 



Z m ~ Z rms = \jk B T jk c . (16) 

As an example, in the experiment 1], for T = 3K, the minimum value of z m is 1.72 x 10 3 . 
Thus, the z-component of the effective field, 2t]z m = 10. Taking e = 2r]z m , we still satisfy 
the adiabatic condition (14) but clearly violate the condition (15) for full reversals. In this 
case, Alj c /u) c = 7 x 10~ 7 , and is 300 times greater than for the experimental parameters 
(13). 

The OSCAR signal should be compared with the frequency thermal noise 5uo c /uo c . To 
estimate this noise we assume that the noise force changes from zero to its rms value F rms 
while the ^-coordinate of the cantilever tip changes from zero to Z m . Then the effective 
change of the cantilever spring constant is 5k c = F rms /Z m . Using the expression uj 2 = k c /m*, 
where m* is the effective cantilever mass, we obtain 

5uj r 5k r F r , 



rms 



u) c 2k c 2k c Z m 
Finally, using the well-known estimate for F rms (see, for example, [6j 



(17) 



F rms = 2(k B Tk c B/Qu c ) 1/2 , (18) 

where B is the measurement bandwidth, we derive the estimate for the frequency thermal 
noise 

' UV {Z m y\k B TB/k c uj c Q) l l 2 . (19) 



(This estimate was derived in |7j, using a different approach.) For the "natural" measure- 
ment bandwidth, B = u c /4:Q, we have from Eq. (19) 

^ = (2Qz m y\k B T/k c y/ 2 . 



7 



Taking parameters (12) and T = 3K, Q = 10 4 , we find 8u c /u c m 9.7 x 10~ 8 , which is smaller 
than the estimated OSCAR signal, Auj c /uj c = 7 x 10~ 7 . 

Note, that the thermal noise is inversely proportional to the cantilever amplitude Z m . 
Thus, the amplification of the OSCAR signal, associated with a decrease of Z m , does not 
assume an increase of the signal-to-noise ratio. However, the decrease of the field amplitude, 
Bi, which becomes possible for small values of Z m , would result in reducing the cantilever 
heating, which is an important obstacle for a single-spin detection. 

So far, we did not discuss the quantum effects in OSCAR. These effects appear if the 
initial average spin is not parallel to the effective field. In the classical case, the angle 9 
between the effective field and the magnetic moment is an integral of motion. The estimated 
classical frequency shift 

ry 1 cos 9 

Au c /u c = - ' (20) 

V 2? 7 Z ra + £ 2 

will approach zero as 9 approaches tt/2. In the quantum case, we expect to obtain one of 
two possible values in (11) independent of the initial angle 9. This outcome corresponds to 
the standard Stern-Gerlach effect. 



IV. SIMULATIONS OF QUANTUM DYNAMICS 

In our simulations we used the following values of parameters 

v = 0.3, z m = 13, e = 10. (21) 

For the cantilever used in Jj| this value of rj corresponds to a magnetic field gradient of 
1.4 x 10 7 T/m, and the value of z m corresponds to the temperature below 170/jlK. 

These parameters allow us to reduce the computational time. The maximum value of 
the ^-component of the effective field is, 2r]z m = 7.8. Thus, the condition (14) of adiabatic 
motion is satisfied while the condition (15) for full reversals is clearly violated. The estimated 
frequency shift is: \Au c \/u c = 7.9 x 10~ 3 . 

The initial conditions which we used correspond to the quasiclassical state of the can- 
tilever, and the electron spin pointing in the positive z-direction. That means, in Eqs. (4), 
up(z, 0) = and, 

oo 

u a (z, 0) = £ n^2 n / 2 a n H n (z) exp[-(z 2 + |« | 2 )/2], (22) 

n=0 



s 




FIG. 2: Functions \u a (z, r)| 2 (black curve) and \u^(z,t)\ 2 (red curve), for six values of r. Initial 
conditions: (z(0)) = z m = 13, {p(0)) = 0, and the electron spin points in the positive z-direction. 



where H n (z) is a Hermitian polynomial, and 

a = ((z(0)) + z(p(0)»/V2, (z(0)) = 13, (p(0)) = 0. (23) 
Our simulations show that the probability distribution for the cantilever position 

P(z,r) = \u a (z,T)\ 2 + \u t3 (z,T)\ 2 (24) 

eventually splits into two peaks which describe two trajectories of the cantilever. The ratio 
of the integrated probabilities for two peaks is given approximately by tan 2 (8/2), where 9 
is the initial angle between the directions of the average spin and the effective field -B e // = 
(e, 0, — 2r)z(0)). The approximate position of the center of the first peak is 

zi = z m cos(l - \ Au c \/u c )t, (25) 

and the approximate position of the second peak is 

z 2 = 2 m cos(l + \Alu c \/u c )t, (26) 
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FIG. 3: Classical and quantum cantilever dynamics in OSCAR MRFM; (a) classical cantilever 
coordinate z(t), (b) quantum average cantilever coordinate (z(t)), (c) the Fourier spectrum for 
(a): z(t) = J2 n A n cos(uj n T + ip n ) , (d) the Fourier spectrum for (b): (z(t)) = J2 n B n cos(u; n T + 6 n ). 
All parameters are the same as in Fig. 2. In a) and b) we show, for convenience, time sequences 
shorter than those used in order to get the Fourier spectrum shown in c) and d). 



where \Auj c \/uj, 



8.0 x 10 



-3 



c ~ u.u ~ xw . (The value \Auj c \/uj c estimated from (11) is almost the 
same, 7.9 x 10~ 3 .) Both functions u a (z,r) and up(z,r) contribute to each peak. (See 
Fig. 2.) When two peaks are clearly separated, the wave function ^ can be approximately 
represented as a sum of two functions ty\ and ^2, which correspond to two peaks in the 
probability distribution. We have found that each function ^1 and ^2, with the accuracy 
to 1%, can be represented as a product of the coordinate and spin functions 



$1 = Ri(z,r) X i(r), ^ 2 = R 2 (z,t) X 2(t). 



(27) 



The first spin function Xi ( T ) describes the average spin which points approximately oppo- 
site to the direction of the effective field -B e //,i = {z, 0, z\), for the first cantilever trajectory. 
The second spin function X2{ T ) describes the average spin which points approximately in 
the direction of the effective field -B e //,2 = (£, 0, z 2 ), for t 

Unlike a conventional MRFM dynamics studied in (3, 
different effective fields for two cantilever trajectories. 

That is why, in general, the average spins corresponding to two cantilever trajectories 



le second cantilever trajectory. 
|^ the OSCAR technique implies 
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FIG. 4: Evolution of the density matrix without the decoherence effects {Q~ l = D = 0), for four 
instants of time (as indicated in the figure) . Parameters: e = 10, rj = 0.3. The left column: 
hi I X)s /°s,sl> the right column: In | J2 S Ps,-s\- Initial conditions: (z(0)) = —8, (p(0)) = 0, and the 
electron spin points in the positive z-direction. Contour lines have been obtained by intersections of 
the functions In | J2s Ps,±s\ with horizontal planes at different heights h. Different colours indicates 
different regions in the following way : h < —14 (white) ,—14 < h < —12 (black) , —12 < h < —10 
(red), -10 < h< -8 (green) -8 < h < -6 (blue) , -6 < h < -4 (yellow), h > -4 (pink) . 

do not point in the opposite directions, and the wave functions Xi( r ) an d X2{ T ) are not 
orthogonal to each other. The only exceptions are the instants r for which two cantilever 
trajectories intersect providing a unique direction for the effective field. 

As we already mentioned, the frequency shift for two cantilever trajectories Au c /u c was 
found to be ±8 x 10~ 3 . For comparison, we solved the classical equations of motion with the 
same parameters and initial conditions as in the quantum case. The frequency shift for a 
single cantilever trajectory in the classical case was found to be Au> c /u c = — 5 x 10~ 3 , which 
is smaller than the quantum shift and very close to the value —5.3 x 10~ 3 derived from the 
estimate (20). (See Fig. 3.) 

V. DECOHERENCE OF THE SCHRODINGER CAT STATE 

To describe qualitatively the decoherence process, we solved numerically the master equa- 
tion (5). The initial density matrix describes a pure quantum state with the quasiclassical 
cantilever and the spin oriented in the positive z-direction. At r = we have 



where the expression for u a (z, 0) is given in (22), and we put in (23) (p(0)) = 0, (-2(0)) = —8. 
Fig. 4 shows the evolution of the density matrix without effects of decoherence, for the 
following values of parameters 



The left column demonstrates the behavior of the spin diagonal density matrix compo- 
nents (contour lines for In | J2 S Ps,s\), an d the right column demonstrates the behavior of the 




(28) 



e = 10, r] = 0.3, Q 



-i 



D 



0. 
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FIG. 5: Evolution of the density matrix: effects of decoherence and thermal noise; D = 20, 
Q^ 1 = 0.001, e = 10, 7] = 0.3. The left column: ln|X) s Ps,s|, the right column: ln|^ s Ps,-s|- 
Initial conditions: (z(0)) = —8, (p(0)) = 0, and the electron spin points in the positive z-direction. 
Contour lines have been obtained by intersections of the functions In | J2 S Ps,±s\ with horizontal 
planes at different heights h. Different colours indicates different regions in the following way : 
h < -14 (white) -14 < h < -12 (black) , -12 < h < -10 (red), -10 < h < -8 (green) 
-8 < h< -6 (blue) , -6 < h < -4 (yellow), h > -4 (pink). 

spin non-diagonal density matrix components (contour lines for In | Y, s Ps,-s\)- Initially we 
have one peak on the plane z — z'. Eventually, this peak splits into four peaks. Two spatial 
diagonal peaks, which are centered on the line z = z', correspond to two cantilever trajecto- 
ries. Two spatial non-diagonal peaks describe a "coherence" between the two trajectories, 
which is a quantitative characteristic of the Schrodinger cat state. All four spin components 
of the density matrix p StS >(z, z' ,t) (with s = ±1/2, s' = ±1/2) contribute to each peak in 
the z — z' plane. 

The density matrix p s , s '(z, z', r) can be represented as a sum of four terms 

p s Az,z',r) = j2pi%,z',T), (29) 
k=i 

where each term describes one peak in the z — z' plane. Suppose that first two terms in 
(29) with k = 1,2 describe the spatial diagonal peaks, and two other terms with k = 3,4 
describe the spatial non-diagonal peaks. 

We have found that the diagonal terms and p^ can be approximately decomposed 
into the tensor product of the coordinate and spin parts 

pW(z, z', r) = RW(z, z', r)xS(r), k = 1,2. (30) 

The spin matrix x^I>{t) describes the average spin which points approximately opposite to 
the direction of the effective field -B e //,i- The spin matrix x s i'( T ) describes the average spin 
which points approximately in the direction of the effective field -B e //,2- 

We have found that the same properties of the density matrix remain for the case 
<5 _1 ,-D 7^ 0. Fig. 5 demonstrates the effects of decoherence and thermal noise for 
Q^ 1 = 0.001 and D = 20. The spatial non-diagonal peaks in the z — z' plane quickly 
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decay. This reflects the effect of decoherence: the statistical mixture of two possible can- 
tilever trajectories replaces the Schrddinger cat state. Next, the spatial diagonal peaks 
spread out along the line z — z' . This reflects the classical effect of the thermal diffusion. 

Conclusion 

We presented a quantum theory of the OSCAR MRFM technique. We demonstrated that 
the OSCAR signal can be significantly amplified by using partial reversals of the effective 
field instead of the full reversals. If the initial angle 9 between the directions of the effective 
field and the average spin is not small, the quantum frequency shift remains constant, 
while the corresponding classical frequency shift is proportional to cos#. This result is a 
manifestation of the Stern-Gerlach effect in the OSCAR MRFM. 

Unlike the "conventional MRFM" in which the reversals of the effective field are pro- 
vided by the frequency modulation of the rf field, the OSCAR MRFM exhibits two possible 
directions of the effective field. As a result, unlike the conventional MRFM, and the Stern- 
Gerlach effect, the two directions of the spin corresponding to two cantilever trajectories do 
not remain antiparallel to each other during the quantum evolution. 
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